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1 Introduction 



The computation of ionization rates or probabilities of atoms through low in- 
tensity (/ << 3.5 x 10 16 W/cm 2 ) laser radiation can be carried out successfully 
using perturbation theory around the solution of the Schrodinger equation with- 
out the presence of the laser fields [|TJ. With the advance of laser technology, 
nowadays intensities of up to (10 19 W/cm 2 ) are possible and pulses may be re- 
duced to a duration of (r ~ 10~ 15 s), [] the region of validity of the above method 
is left. The new regime is usually tackled by perturbative methods around the 
Gordon- Volkov solution || of the Schrodinger equation H, |[ |7], |], |9|, |10| , fully 
numerical solutions of the Schrodinger equation JTT], [12], [13], [14], [TJ|, [RJ , Floquet 



solution |17], [18|, [L9| , high frequency approximations p0[ or analogies to classical 
dynamical systems All these methods have its drawbacks. The most sur- 
prising outcome of the analysis of the high-intensity region for short pulses (the 
pulse length is smaller than 1 ps) is the finding by the majority of the atomic- 



physics community (see [22], |23], |24|, |25| and references therein) of so-called atomic 
stabilization. This means that the probability of ionization by a pulse of laser 
radiation, which for low intensities increases with increasing intensities reaches 
some sort of maximum at high intensities and commences to decrease until ion- 
ization is almost totally suppressed. This picture is very counterintuitive and 
doubts on the existence of this phenomenon have been raised by several authors 
|9|, |26], [27], |7|, |1(J, who do not find evidence for it in their computations. So far no 
support is given to either side by experimentalists.^] For reviews on the subject 



we refer to 22, 23, 24, 25 



Since all of the above methods involve a high degree of numerical analysis, 



*For a review and the experimental realization of such pulses see for instance [|]. 

^ Experimental evidence for some sort of stabilization is given in pqj , but these experiments 

deal with intensities of 10 13 W/cm 2 , which is not the " ultra- intense" regime for which the 

theoretical predictions are made. 
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which are difficult to be verified by third parties, it is extremely desirable to 



reach some form of analytical understanding. In |29], |30|, [U], |32| we derived an- 
alytical expressions for upper and lower bounds for the ionization probability, 
meaning that the ionization probability is certainly lower or higher, respectively, 
than these values. The lower bound in particular may be employed to investigate 
the possibility of stabilization for an atomic bound state. In [[32] we analyzed 
the hydrogen atom and found that for increasing intensities the lower bound also 
increases and hence that the existence of atomic stabilization can be excluded in 
the sense that the ionization probability tends to one. The shortcoming of our 
previous analysis []32| is, that definite conclusions concerning the above question 



may only be reached for extremely short pulses (r < 1 a.u.), which are exper- 
imentally unrealistic. In the present article we analyze these bounds in further 
detail and demonstrate that atomic stabilization can also be excluded for longer 
pulses. 



Some authors || [14]] put forward the claim that in order to "observe" atomic 
stabilization one requires pulses which are switched on, sometimes also off, smoothly. 
This seems very surprising since stabilization is supposed to be a phenomenon 
specific to high intensities and with these type of pulses emphasis is just put on 
the importance of the low intensity regime. It further appears that among the 
authors who put forward these claims, it is not commonly agreed upon, whether 
one should associate these pulse shapes to the laser field or to the associated 
vector potential. We did not find a proper and convincing physical explanation 
why such pulses should produce so surprising effects in the literature. Geltman 
[T0|| and also Chen and Bernstein [2(| do not find evidence for stabilization for 
these type of pulses with smooth and turn on (and off) of the laser field. 

In order to address also the validity of these claims in our framework we extend 
in the present paper our previous analysis to various type of pulses commonly 
employed in the literature in this context and investigate also the effects different 
frequencies might have. Once more we conclude that our arguments do not 



support atomic stabilization. 

Our manuscript is organized as follows: In section 2. we briefly recall the 
principle of our argumentation and our previous expressions for the upper and 
lower bounds for the ionization probability and discuss them in more detail for the 
hydrogen atom. We then turn to an analysis for specific pulses. In section 3. we 
state our conclusions. In the appendix we present the explicit computation for the 
Hilbert space norm of the difference of the potential in the Kramers-Henneberger 
frame and the one in the laboratory frame for any bound state. 



2 The upper and lower bounds 

For the convenience of the reader we commence by summarizing briefly the main 
principle of our argument. Instead of calculating exact ionization probabilities 
we compute upper and lower bounds for them, meaning that the exact values are 
always greater or smaller, respectively. We then vary these bounds with respect 
to the intensity of the laser field and study their behaviour. If the lower bound 
tends to one with increasing intensity, we can infer that stabilization is definitely 
excluded. On the other hand, if the upper bound tends to zero for increasing 
intensities, we would conclude that stabilization is present. In case the lower 
bound increases, but remains below one, we only take this as an indication for 
a general type of behaviour and interpret it as not providing any evidence for 
stabilization, but we can not definitely exclude its existence. In case the lower 
bound becomes negative or the upper bound greater than one, our expressions 
obviously do not allow any conclusion. 

The non-relativistic quantum mechanical description of a system with po- 
tential V in the presence of linearly polarized laser radiation is given by the 
Schrodinger equation involving the Stark Hamiltonian 

.dip(x,t) ( A 

i - 



dt 



= (-- + V + z-E(t)^(x,t) = H(t)^(x,t). (2.1) 
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For high, but not relativistic, intensities the laser field may be approximated 
classically. We furthermore assume the dipole approximation. In the following 
we will always use atomic units h = e = m e = c ■ a = 1. For a general time 
dependent Hamiltonian H(t) the ionization probability of a normalized bound 
state ip is defined || |29[ as 



P(V0 = ||(1-P + )^|| 2 = 1-||7V^H 2 • (2-2) 

The gauge invariance of this expression was discussed in [[52]. Here \\if>\\ denotes 
as usual the Hilbert space norm, i.e. HV'II 2 = (ip,^) = I \ip{x)\ 2 d 3 x. We always 
assume that H± = lim H(t) exists and ip is then understood to be a bound 

t— >±oo 

state of H-. V + and V- denote the projectors onto the space spanned by the 
bound states of H + and H- , respectively and S is the unitary "scattering matrix" 

S= lim exp(it + H + ) -C/(t+,t_) ■ exp{-it_H_) . (2.3) 

Here the unitary time evolution operator U(t+,t-) for H(t), brings a state from 
time t- to t + . Note that by definition < P(ip) < 1. Employing methods of 



functional analysis we derived in |29], [30], |31|, ^2| several analytical expressions by 
which the possible values for the ionization probability may be restricted. We 
emphasize once more that these expressions are not to be confused with exact 
computations of ionization probabilities. We recall here the formula for the upper 



PuW*= J \\(V(x-c(t)e z )-V(x))^\\dt+\c(r)\ \\p z ^\\ + \b(r)\ \\z^\\ (2.4) 
o 

and the lower bound 

PiW = l - [j \\{V{x-c{t)e z )-V{x))^\\dt (2.5) 

2 , , . , 2|6(r)| 



\\(V(x-c(r)e z ) - V(x)W\ + n J\ \J Pz 



2E + b(r) 2U " v v ' " v JJTn 2£ + 6(r) 2 ' } 
which were deduced in J3"2"| . e z is the unit vector in the z-direction. Here we use 
the notation 

Ut) := f E(s)ds c(t) := /* b(s)ds , (2.6) 

Jo Jo 



for the total classical momentum transfer and the total classical displacement, 
respectively. Note that for the vector potential in the z-direction we have A(t) = 
— \b{t) + const . It is important to recall that the expression for the lower bound is 
only valid if the classical energy transfer is larger than the ionization energy of the 
bound state, i.e. \b 2 {r) > —E. Our bounds hold for all Kato small potentials!] 
In particular the Coulomb potential and its modifications, which are very often 
employed in numerical computations, such as smoothed or screened Coulomb 
potentials, are Kato small. However, the delta-potential, which is widely used in 
toy-model computations because of its nice property to possess only one bound 
state, is not a Kato potential. 

In the following we will consider a realistic example and take the potential 
V to be the Coulomb potential and concentrate our discussion on the hydrogen 
atom. In this case it is well known that the binding energy for a state ip n im is 
E n = -2^2, IbWViOoll 2 = 3^2 and ||^noo|| 2 = |(VV)ok 2 |VV)o) = if(5n 2 + 1) (see 
for instance |33]). We will employ these relations below. In [[32| it was shown, 



that the Hilbert space norm of the difference of the potential in the Kramers- 
Henneberger frame pi [37] and in the laboratory frame applied to the state if> 



N(y,if;) := \\{V{x-y)-V{x))i,\\ (2.7) 

is bounded by 2 when if) = if) 1Q0 for arbitrary y = ce z . We shall now investigate 
in more detail how this function depends on c. In order to simplify notations we 
ignore in the following the explicit mentioning of e z . In the appendix we present 
a detailed computation, where we obtain 

N 2 (c,ip W0 ) = 2 + (l + |c|- 1 )e-l c l^(|c|) + (l-| C |- 1 )el c l J Ez(-|c|) + ^(e- 2 l c l - l) . 

' (2-8) 

* Potentials are called Kato small if for arbitrary there < a < 1 there is a constant b < oo, 
such that \\Vtp\\ < a\\ — Aip\\ + b\\ip\\ holds for all tp in the domain T>(Hq) of Hq = —A/2, see 
for instance [B4L M . 
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Here Ei(x) denotes the exponential integral function, given by the principal value 
of the integral 



oo . 

e 1 



Ei(x) = - J — dt for x > . (2.9) 

— X 

Considering now the asymptotic of N, we obtain as expected limiV = and 

c^O 

lim N = y/2. Noting further that iV is a monotonically increasing function of c, 
(one may easily compute its derivatives w.r.t. c, but we refer here only to the 



plot of this function in figure 1), it follows that our previous |32j estimate may 
in fact be improved to N(c, "0ioo) < v2. The important thing to notice is, that 
since N(c, ipioo) is an overall increasing function of c, it therefore also increases 
as a function of the field strength. The last term in the bracket of the lower 
bound Pi(ip) is a decreasing function of the field strength, while the second term 
does not have an obvious behaviour. Hence if the first term dominates the whole 
expression in the bracket, thus leading to a decrease of Pi{ip), one has in principle 
the possibility of stabilization. We now investigate several pulse shapes for the 
possibility of such a behaviour and analyze the expressions 



P„(<>„i = >^ I N(c(t),ip mi )dt+'-^ + \b(T)\} . (2.II) 




Here we have simply inserted the explicit values for E\, ll^iooll and \\p z ip 



100 



into (|2.4|) and (|2.5|) , and understand ^(0,^100) to be given by the analytical 
expression (|2.8| ). The formulae presented in the appendix allow in principle also 
the computation of N(c, ipnim) for different values of n, I and m. However, for 
I 7^ the sum over the Clebsch-Gordan coefficients becomes more complicated 
and due to the presence of the Laguerre polynomial of degree n in the radial 
wave-function R n \ this becomes a rather complex analytical computation. We 
will therefore be content with a weaker analytical estimate here. In fact, we have 



iV 2 (c(t),Vnoo) < 2(iP n00 ,V(x) Vnoo) = ^ • (2.12) 



In the appendix of 3J3] this statement was proven for n — 1. The general proof 
for arbitrary n may be carried out exactly along the same line. Therefore, we 
obtain the following new upper and lower bounds 

[ 2 4 1 1 2|6(r)| I 2 , 

PUM - 1 - j-^r + 6(r)2 _ — + ^ 6(r)2 _ 1/n2 j (2-13) 



2 |c(t)| /5n 2 + 1 ' 2 



*WVW») = |^ r+ ^7f + O - |6(r)l | ' (2 ' 14) 



which are weaker than Q2.11D and (|2.10|) , in the the sense that 



Plw(^nOo) < Pl{lpn00) < -P(^nOo) < -P«(^nOo) < Puw^nOo) ■ (2.15) 

In order for fl2.13| ) to be valid we now have to have b(r) 2 > \. We will now turn 
to a detailed analysis of these bounds by looking at different pulses. Our main 
purpose in the present manuscript for considering states of the type ip n i m with 
n 7^ is to extend our discussion to pulses with longer duration, see also section 
2.3. The reason that longer pulse durations are accesible for states with higher n 
is the n-dependence in estimate (|2.15| ) and its effect in fl2.14|) and ( |2.13j ). 



2.1 Static Field 

This is the simplest case, but still instructive to investigate since it already con- 
tains the general feature which we will observe for more complicated pulses. It 
is furthermore important to study, because it may be viewed as the background 
which is present in most experimental setups, before more complicated pulses can 
be generated. For a static field of intensity / = Eq we trivially have 

E t 2 

E{t) = E b{t)=E t c (t) = _£- (2.16) 

for < t < t. Inserting these functions into ( |2.10p we may easily compute the 
upper and lower bound. Here the one dimensional integrals over time, appearing 
in ( |2.11 ) and ( p,10| ) were carried out numerically. The result is presented in figure 



2, which shows that a bound for higher intensities always corresponds to a higher 



ionization probability. The overall qualitative behaviour clearly indicates that 
for increasing field strength the ionization probability also increases and tends to 
one. In particular lines for different intensities never cross each other. Surely the 
shown pulse lengths are too short to be realistic and we will indicate below how to 
obtain situations in which conclusive statements may be drawn concerning longer 
pulse durations. In the following we will always encounter the same qualitative 
behaviour. 

2.2 Linearly polarized monochromatic light (LPML) 

Now we have 



for < t < t. The result of the computation which employs these functions 



bounds indicate that for increasing field strength the ionization probability also 
increases. Keeping the field strength fixed at E = 2 a.u., a comparison between 
the case for uj = 0.4 and to = 4 shows (figure 4), as expected, the lower bounds for 
the ionization probability to be decreasing functions of the frequency. The peak 
on the left, which seems to contradict this statement for that region, is only due 
to the fact that the expression for the lower bound is not valid for uj = 0.4 in that 
regime. Clearly, this is not meant by stabilization, since for this to happen we 
require fixed frequencies and we have to analyze the behaviour for varying field 
strength. The claim |TJ], |2D| is that in general very high frequencies are required 
for this phenomenon to emerge. Our analysis does not support stabilization 
for any frequency. As mentioned above, the shortcoming of the analysis of the 
bounds -P«(V'ioo) an d -Pzt^ioo) is that we only see an effect for times smaller than 
one atomic unit, figure 4 and figure 5 also show that by considering -P(Vvoo) 
for higher values of n our expressions allow also conclusions for longer pulse 






(2.17) 



in order to compute ( |2.1U| ) and ( |2.11| ) is illustrated in figure 3. Once again our 
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durations. For the reasons mentioned above, in this analysis we employed the 
slightly weaker bounds Q2.14J) and Q2.13|) . 



2.3 LPML with a trapezoidal enveloping function 

We now turn to the simplest case of a pulse which is adiabatically switched on and 



off. These type of pulses are of special interest since many authors claim [14], || 
that stabilization only occurs in these cases. We consider a pulse of duration t$ 
which has linear turn-on and turn-off ramps of length T. Then 



E(t) = E sm(ut) 



± for < t < T 



1 for T<t<(r -T) (2.18) 
^ for (r - T) < t < t 

6(r ) = 4^{ sin (^ T )- sin (^ r o) + sin (^( r o- T ))} (2-19) 



E 



( 



c(to) = —5 — 2 — 2 cos(c<jT) + 2 cos(co>t ) — 2 cos(u;(t — T)) 
uj a T \ 

—uTsin(ujT)+UTosm(uT)+uTsin(oj(ro — T)^ . (2.20) 

The expressions for b(t) and c(t) are rather messy and will not be reported here 
since we only analyze the weaker bounds. Notice that now, in contrast to the 
previous cases, both 6(t ) and c(t ) may become zero for certain pulse durations 
and ramps. We shall comment on this situation in section 3. We choose the 
ramps to be of the form T = (m + ^) 77 ( m being an integer) for the lower and 
T = (m + — for the upper bound. Our lower bound does not permit the 



2, 

analysis of half cycles since then b(r ) = 0. The results are shown in figure 6 and 
7, which both do not show any evidence for stabilization. They further indicate 
that a decrease in the slopes of the ramps with fixed pulse duration, leads to a 
smaller ionization probability. Once more (we do not present a figure for this, 
since one may also see this from the analytical expressions), an increase in the 
frequency leads to a decrease in the lower bound of the ionization probability for 
fixed field strength. 
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2.4 LPML with a sine-squared enveloping function 

Here we consider 



E(t) = E sin 2 (ttt) sin(ut) (2.21) 

b(t) = [ 8 n 2 + 2 uj 2 cos(c; t) - 8 Q 2 cos(cj t) 

16 uj ii — 4ur \ 

-uj 2 cos((cj -2Q) t) -2uj Q cos((u; -2Q)t) 

-uj 2 cos((cj + 2 Q) t) + 2 Q cos((cj + 2 J?) t)J (2.22) 

c(t) = 5, k \ -8oj 3 n 2 t + 32un 4 t-2u 4 sm(ut) 

W 4cu 2 (cu-2/2) 2 (cj + 2^) 2 V 

+16u; 2 J? 2 sin(u;i) — 32 J? 4 sin(u;t) — co> 4 sin((2 Q — u) t) 

-4uj 3 Q sin((2 - uj) t) - 4w 2 Q 2 sin((2 - w) t) 

+uj 4 sin((u; + 2 fl) t) - 4cj 3 /2 sin((u; + 21?) t) 

+4 aj 2 Q 2 sin((cj + 2/2) t) J (2.23) 

for < t < t. At first sight it appears that both b(t) and c(t) are singular at 
uj = ±2Q, which of course is not the case since both functions are bounded as one 
may easily derive. With the help of the Schwarz inequality it follows that always 



\b(t)\ < t^ll^H and \c(t)\ < ~ti||.E||. We first investigate the situation in which 
this pulse is switched on smoothly but turned off abruptly. Figure 8 shows that 
the bounds become nontrivial for times larger than one atomic unit in the same 
fashion as in the previous cases by considering Pi(i/j n oo) for higher values of n. 
Figure 9 shows that also in this case the ionization probability tends to one and 
no sign for stabilization is found. Figure 10 shows the lower bound in which the 
pulse length is taken to be a half cycle of the enveloping function. Once more it 
indicates increasing ionization probability with increasing field strength and also 
for increasing values for n. Following now Geltman |10| and Su et al. [1JJ] we 



employ the sine-square only for the turn-on and off and include a plateau region 
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into the pulse shape. Then 



E(t) = E sm(ujt) 



Kro) 
c ( r o) 



sin' 
1 

sin' 



for < t < T 
for T <t< (t - T) 
(ZEM) for (T -T)<t<T 



E 7t 2 (1 + cos(u T) — cos(c<j (T — To)) — cos(u T )) 



Enn 2 2u 2 



{7T . - W «T 2 )' 

_2 



2ujti 2 -2uj 3 T 2 

_2 _ , ,3 t-i2 



(2.24) 



(2.25) 



I ^' 7r J To — uj 3 T 2 To — ujtt 2 T cos(u; T) + uj 3 T 3 cos(c<j T) 
+ uj n 2 t cos(u; T) — u; 3 T 2 r cos(u; T) — u it 2 T cos(u; (T — t )) 



+ uj 3 T 3 cos(cj (T - r )) + 7T 2 sin(a; T) - 3 uj 2 T 2 sin(cj T) 

+ 7r 2 sin(u; (T — t )) — 3u 2 T 2 sin(u; (T — t )) — n 2 sm(uT ) 

+ 3cu 2 T 2 sm(LUT )) . (2.26) 

(Also in these cases the apparent poles in 6(r ) and c(r ) for uj = ±^ are ac- 
companied by zeros.) The results of this computations are shown in figure 6 and 
7, once more with no evidence for bound-state stabilization. A comparison with 
the linear switch on and off shows that the ionization probability for sine-squared 
turn-on and offs is lower. The effect is larger for longer ramps. 



3 Conclusions 

We have investigated the ionization probability for the hydrogen atom when 
exposed to ultra-intense shortly pulsed laser radiation of various types of pulse 
shapes. In comparison with |32| , we extended our analysis to the situation which 



is applicable to any bound-state ip n i m and in particular for the Vioo-state we 
carried out the computation until the end for the stronger upper ([2.4]) and lower 
bounds. We overcome the shortcoming of [32] which did not allow definite 



statements for pulses of durations longer than one atomic unit by investigating 
the bounds for higher values of n. A direct comparison between existing numerical 
computations for small n, in particular n=l, and reasonably long pulse durations 
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is at present not feasible. As our computations show (see also J331) there is 



of course a quantitative different behaviour for different values of n. However, 
qualitatively we obtain the same behaviour (refer figure 10) and therefore we do 
not think this to be of any physical significance. It would be very interesting to 
carry our analysis further and also investigate the effect resulting from varying 1 
and m. In principle our equations already allow such an analysis, but due to the 
sum in (A. 6) the explicit expressions will be rather messy and we will therefore 
omit them here. 

We regard the lack of support for the existence of bound-state stabilization 
in a realistic three dimensional atom resulting from these type of arguments, 
even for high values of n, as more convincing than for instance the support for 
stabilization based on one-dimensional toy models. 

For the situation when the total classical momentum transfer b{r) and the to- 
tal classical displacement c(r) are non-vanishing we confirm once more the results 
of |32j and do not find any evidence for bound-state stabilization for ultrashort 



pulses. This holds for various types of pulses, whether they are switched on (and 
off), smoothly or not. We therefore agree with Geltman in the conclusion that 
smooth pulses in general will only prolong the onset of ionization but will not 
provide a mechanism for stabilization. 

There is however a particular way of switching on and off, such that b{r) = 0, 
but c(r) 7^ 0. These type of pulses are used for instance in || [14]]. Unfortunately 
our bounds do not permit to make any definite statement about this case, since 
the lower bound is not applicable (in the sense that then the necessay condition 
\b 2 (r) > —E for th validity of the lower bound is not fulfilled) and the upper 
bound gives for typical values of the frequency and field strength ionization prob- 
abilities larger than one. So in principle for these type of pulses the possibility of 
bound-state stabilization remains. It would be very interesting to find alternative 
expressions for the upper and lower bound which allow conclusions on this case. 

For the case 6(r) = c(r) = the upper bound P u remains an increasing 
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function of the field strength due to the properties of the Hilbert space norm 
of the difference of the potential in the Kramers- Henneberger frame and in the 
laboratory frame applied to the state V'ioo- The weaker upper bound takes on the 
value Puwi^noo) = ~f) which at first sight seems counterintuitive, since it implies 
that the upper bound decreases with increasing n, i.e. for states close to the 
ionization threshold, and fixed r. Classically this may, however, be understood 
easily. For closed Kepler orbits, i.e. ellipses, with energies sufficiently close to 
zero (depending on r), for any pulse with small b(r) and c(r), these quantities 
will be very close to the actual changes, caused by the pulse, of the momentum 
and the coordinate, respectively. So in this case ionization, i.e. the transition 
to a hyperbolic or parabolic orbit will therefore be very unlikely. This kind of 



behaviour was also observed in |38| for a Gaufiian pulse, for which 6(r) = and 
c(r) + 0. 

Acknowledgment: We would like to thank J.H. Eberly, S. Geltman and V. 
Kostrykin for very useful discussions and correspondences. One of the authors 
(CFMF) is supported by the DAAD. 

Appendix 

In this appendix we will provide the explicit calculation of the term 

N 2 (y,^j) = (4>,V(x) 2 ij) + (^,V(x-y) 2 ^)-2(4>,V(x-y)V(x)^) (A.l) 



For ip = ipnim the first term is well known to equal ^wn^y^ p3| . We did not find 
a computation for the matrix element involving the Coulomb potential in the 
Kramer s-Henneberger frame in the literature and will therefore present it here. 
Starting with the familiar expansion of the shifted Coulomb potential in terms 
of spherical harmonics 



' t(4l)\h^7^o^A) (A.2) 



\z-y\ fa\ r > V V 2Z + 1 
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where r< = Min(\x\, \y\) and r> = Max(\x\, \y\), we obtain 

°° r r 

(*nlm\ \x-y\- 1 Ifl-^n/m) = £ / dAY,*^ 0*5 

l'=0 



47T 



2/' + 1 



which by the well known formula from angular momentum theory 

jny* ~y ~y 

I m 1 h mi 1 h r?i2 



\ {2h 4n(2l+\) ^ {lll2] °° m {kh] miTn2 ^ (A ' 3) 
leads to 

OO / 1^ / \ '' + 1 °° / I -H \ \ 

£(« / ;00|ZO)(M , ;mO|Zm> (/ rfr (jt) + / dr i^f j R 2 nl J . (A.4) 

Here (hh] mim2\lm) denote the Wigner or Clebsch-Gordan coefficients in the 
usual conventions (see e.g. |39|] ). 
We shall now consider the term 

(^nlm\ \x-y\~ 2 \^>nlm) (A.5) 

Employing (|A.2|) and the formula 



Y hmi Y hm2 = y (2/l + l) ^ h + (^pi)^^' m ^\rm')(hi 2] ooiro) 



yields 



E E "TO J^^ °l /0 ) 2 no • (A.6) 



l--^ 2 Uut-^> +v+2 V2/ + 1 



Once again applying ( |A.3|) shows that ( |A.5| ) equals 

I' +1+2 00 , M > i'+Z 



00|[0) 2 (r/;0m|/m)(/l;00|/0) ( j dr R 2 nl + J dr R 2 nl 

(A.7) 
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For s-states, i.e. (1 = 0), we may carry out the sums over the Clebsch-Gordan 
coefficients easily. In ( [A .4 ) the only contribution comes from V = and we 
trivially obtain 

(# n00 | \x-y\~ 1 Ifr^noo) = J dr-^R 2 n0 + J drR 2 n0 . (A.8) 



In flA~?l ) the sum over I contributes only for I = and together with (IV; 00 1 00) " 



JS^- it leads to 



(^ n00 | \x- 



21+2 
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^] Rl l+ ldr[^) R 



(A.9) 



We turn to the case n = 1 (with \I/ioo = f° r which ( |A.4j ) becomes 



(^ioo| \x~y\ 1 |*ioo) 



1 - e' 2 ^ 



(A.10) 



As consistency check one may consider the asymptotic behaviours \y\ ^ oo and 
— > 0, which give, as expected, and 2 respectively. Using the series expansion 
for the logarithm, ( |A.9| ) for n = 1 becomes 



(^ioo I |aH7| 2 |*100> = 

Using then the integrals 
J dr ln(l ± r)re 



\y I 



c?r In 



1 2/ 1 + r 



re 2r + I dr In 



r+Jy| 



re 



-2r 



(A.H) 



2cr 



— (jl T 2c)e ±2c ^( T 2c(l±r)) 
-e- 2cr (! + (! + 2cr) ln(l ± r)) j 



(A.12) 



y dr ln(l ± r 1 )re 



l\„„-2cr 



4c 2 



(1 T 2c)e ±2cr £z (2c(Tl-r)) 



-Ei(-2cr) - (1 + 2cr)e- 2cr ln(l ± r" 1 )^ (A.13) 
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we obtain 



(tfiool \x-y/2\- 2 |* 100 > = (1 - W^e-^mm + (1 - \y\~ l )e^Ei{-\y |) 

(A. 14) 

As a consistency check we may again consider the asymptotic behaviour, that is 
\y\ — > and \y\ —>■ oo, which gives correctly 2 and 0, respectively. Assembling 
now ( A.l ), ( A.10|) and (|A.14|) gives as claimed ( |2.8|) . In the same fashion one 
may also compute N(y, ip n i m ) for arbitrary n, I and m. 
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Figure 1: The Hilbert space norm of the difference of the potential in the Kramers- 
Henneberger frame and in the laboratory frame applied to the state ^100 versus 
the classical displacement c. 

Figure 2: Upper (three curves on the left) and lower bound (P t and P u ) for the 
ionization probability of the ^ 10 o-state through a static laser field Eq. The dotted 
line corresponds to E = 5 a.u., the dashed line to E = 10 a.u. and the solid 
line to Eq = 20 a.u. The time is in a.u. 

Figure 3: Upper (three curves on the left) and lower bound (Pi and P u ) for the 
ionization probability of the ^loo-state through a linearly polarized monochro- 
matic laser field E(t) = E sm(ut); oj = 1.5 a.u. The dotted line corresponds to 
E = 5 a.u., the dashed line to E = 10 a.u. and the solid line to E = 20 a.u. 
The time is in a.u. 

Figure 4: Lower bound (Pi w ) for the ionization probability of the ^looo-state 
through a linearly polarized monochromatic laser field E(t) = E sin(ujt), E = 2 
a.u. The dotted line corresponds to u — 0.4 a.u. and the solid line to to = 4 a.u. 
The time is in a.u. 

Figure 5: Lower bound for the ionization (Pi w ) probability of the ip2o oo-state 
through a linearly polarized monochromatic laser field E(t) = E sm(ut), ou — 1.5 
a.u., Eq = 20 a.u. The time is in a.u. 

Figure 6: Lower bound (P/ w ) for the ionization probability of the ■^34 o-state 
through a linearly polarized monochromatic laser field with a trapezoidal and 
a sine-squared turn-on and turn-off enveloping function, upper and lower curve 
of the same line type, respectively, (solid line: | — 12 — | pulse, dashed line: 
I — 10 — I pulse and dotted line: ^ — 6 — ^ pulse), u — 1.5 a.u. 
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Figure 7: Upper bound (Pi w ) for the ionization probability of the ^34 oo-state 
through a linearly polarized monochromatic laser field with a trapezoidal and a 
sine-squared turn-on and turn-off enveloping function, upper and lower curve of 
the same line type, respectively, (solid line: | — 6 — \ pulse, dashed line: § — 4 — | 
pulse and dotted line: | — 2 — | pulse), a; = 1.5 a.u. 

Figure 8: Lower bound (Pi w ) for the ionization probability of the ^30 oo-state 
through a linearly polarized monochromatic laser field with a sine-squared en- 
veloping function E(t) = E sin(ujt) sin(fit) 2 , u = 0.2 a.u., Q = 0.01 a.u., E = 20 
a.u. The time is in a.u. 

Figure 9: Lower bound (Pi w ) for the ionization probability of the ip^ oo-st at e 
through a linearly polarized monochromatic laser field with a sine-squared en- 
veloping function E(t) = E sin(u;t) sin(f2t) 2 , uo = 0.2 a.u., Vl = 0.01 a.u. The 
dotted line corresponds to E = 5 a.u., the dashed line to E = 10 a.u. and the 
solid line to E = 20 a.u. The time is in a.u. 

Figure 10: Lower bound (Pi w ) for the ionization probability of the ^ n oo-state 
through a linearly polarized monochromatic laser field with a sine-squared en- 
veloping function E(t) = E sm(ut) sin(Qt) 2 , u = 0.8 a.u., Q = cj/13.5 a.u. The 
dotted line corresponds to n = 40, the dashed line to n = 35 and the solid line 
to n = 30. 
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